Nucleosynthetic Constraints on the Mass of the Heaviest Superaovae 



Justin M. Brown^ and S. E. Woosley^ 



ABSTRACT 

Assuming a Salpeter initial mass function and taking the solar abundances as a 
representative sample, we explore the sensitivity of nucleosynthesis in massive stars 
to the truncation of supernova explosions above a certain mass. It is assumed that stars 
of all masses contribute to nucleosynthesis by their pre-explosive winds, but above a 
certain limiting main sequence mass, the presupemova star becomes a black hole and 
ejects nothing more. The solar abundances from oxygen to atomic mass 90 are fit 
quite well assuming no cut-off at all, i.e., by assuming all stars up to 120 Mq make 
successful supemovae. Little degradation in the fit occurs if the upper limit is reduced 
to 25 M0. The limit can be further reduced, but the required event rate of supemovae 
in the remaining range rises rapidly to compensate for the lost nucleosynthesis of the 
more massive stars. The nucleosynthesis of the 5-process declines precipitously and 
the production of species made in the winds, e.g., carbon, becomes unacceptably large 
compared with elements made in the explosion, e.g., silicon and oxygen. However, by 
varying uncertain physics, especially the mass loss rate for massive stars and the rate 
for the ^^Ne(a,n)^^Mg reaction rate, acceptable nucleosynthesis might still be achieved 
with a cutoff as low as 18 Mq. This would require a supernova frequency three times 
greater than the fiducial value obtained when all stars explode in order to produce the 
required '^O. The nucleosynthesis of ^°Fe and ^^Al is also examined. 

Subject headings: Nucleosynthesis, Massive Stars, Stellar Evolution, Initial Mass 
Function, Galactic Enrichment 



1. Introduction 

Just which massive stars explode as supemovae and which collapse to black holes has been a 
topic of great interest for a long time. In particular, stellar nucleosynthesis can be used to constrain 
the maximum mass of the supernova that needs to (or ca n) explode in order to explain the abun 



dances of the elements seen in the Sun and elsewhere (e.g. lTwarog & Wheeleiill982 
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19921). Previous works have generally focused on the production of helium and oxygen and the 
ratio A y / AZ where Y is the helium mass fraction and Z is the heavy element fraction. iMaeder 
(Il992|) . for example, concludes that the observed abundances are best fit if stars of all masses con- 



tribu te their pre-exp losive winds, but only stars below 20 - 25 explode as supemovae (though 
see Prantzosl (1 1994|) who contests this conclusion). The remainder presumably end as black holes 



which accrete the remaining star, including all its heavy elements. 

Added interest in this issue has been generated recently by the increasin gly tight observational 
constraints placed upon the masses of presupe mova stars. 



Smartt (2009) finds no evidence for 



supernova progenitors with masses over 20 M0. iHoriuchi et al.l (|201lh also finds an inconsistency 
with the measured rate of core collapse supemovae and the cosmic star formation rate in the sense 
that more stars seem to form than are observed to die as supemovae by a factor of about 2. Of 
course, these arguments are not yet absolute. More massive supemovae may be hidden in dust 
and the connection between main sequence mass and presupemova mass relies on theory. Star 
formation rates and supemova rates are not precisely known, but these constraints may become 
tighter with time and certainly suggest that not all massive stars make luminous supemovae. 

On the theoretical front, it has bee n known for a long time that more massive stars are harder 
to blow up than lower mass ones (e.g. lFryerlll999t Fryer & Hegerll2000|) . With higher mass, the 
entropies around the collapsing iron cores are greater, and the fall off of density with radius is 
consequently more gradual. During the collapse this means a greater accretion rate on the proto- 
neutron star that is harder to reverse. [O'Connor & Ottld201l|) recently quantified this effect with a 
"compactness" parameter (see their Fig. 9) and measured the difficulty of blowing up stars in a ID 
code as a function of that parameter. Their conclusion, interestingly again was that stars heavier 
than 20 M0 were harder to explode. 

These considerations motivate a revisit of the problem of stellar nu cleosynthesis as a functio n 
of mass. Using a fiducial model set of solar metall icity models from IWoosley & Hegeii (|2007|) . 
the same model set used by lO' Connor & OttI (|201l|) . we study the nucleosynthesis and remnant 
masses resulting if the supemov a explosions are truncated above a certain mass, Mbh- A Salpeter 
initial mass function is assumed (|Salpeterlll955|) with a slope of F = -1 .35. We determine not only 
the bulk nucleosynthetic properties like oxygen and remnant mass as a function of Mbh, but also 
examine the synthesis of the i'-process, the individual ratios of important intermediate mass and 
light elements, and the synthesis of interstellar radioactivities, ^^Al and ^"Fe. 
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2. Models 



The yield tables of IWoosley & Hegeii (l2007h give the nucleosynthesis of all species from 
hydrogen through lead for supemovae resulting from non-rotating massive stars with the following 
initial masses: 12-33 (every integer mass), 35-60 (every five masses), 60-80 (every 10 masses), 
100, and 120 solar masses. The authors calculated explosions for four sets of models parameterized 
by the mass cut and explosion energy. Here we use their standard set for which the explosion 
energy was 1.2 x 10^' erg and the mass cut was located at the " entropy jump" whe re S/NAk = 
4.0. These are their "A" mode ls and are the sanie mode ls for which IZhang et al.l (|2008|) calculated 
compact remnant masses and lO' Connor & OttI (|2011|) studied compactness. The spread of these 
models can be seen in Fig. [TJ Nucleosynthesis ejected by the pre-explosive win ds and in the 
explosions was archived separately but is not available in IWoosley & Hegeii (|2007|) (see Figure 4 
in their paper for the mass loss prescription used). The values of some key species are provided 
in an online supplement for this paper. Using this grid of nucleo syntheti c yields, we constructe d a 
stellar population using the high-end initial mass function described by iReid & WilsonI (|2006|) as 
illustrated in Fig. [U 



$(M)d(M) oc M"^-^^^"-^ d(M) (1) 

The total yields of the stellar population were calculated by integrating the yields over the mass 
function, as described in eq. dU, where m, is the total production (in solar masses) of isotope /, and 
Ei(M) and Wi(M) are the total ejecta of isotope / from, respectively, the supernova explosion and 
the winds in solar masses from a star with initial mass M: 

rM upper I-UOMq 

mi= / <!>(M)Ei(M)dM+ / $(M)W^(M)dM (2) 

JuMq JuMq 

We calculate the yields for the stable isotopes, varying the upper mass limit to the values listed in 
Table [B 

Our results are expressed in terms of a simple "production factor" for each isotope defined as 

m,/ V . nij 

Pi= , (3) 



where S, is the mass fraction of the isotope in the Sun (iLoddersI I2003h . The importance of the 
production factor lies in its relationship to the supernova rate. It can be shown that the production 
factor from high-mass progenitors is inversely proportional to the number of supemovae estimated 
by the model, assuming that the majority of the isotope's production comes from this stellar popu- 
lation. 



For comparison, we plot the produ ction factors as a functio n of atomic mass for a population 
with no upper bound in Fig. |2](see also lWoosley & Hegeiil2007|) . 
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3. Limits on Mbh 

The mass of the heaviest supernova that has to explode is constrained by a variety of obser- 
vations. These include not only the nucleo synthetic pattern of the intermediate mass elements, but 
especially the light component of the i'-process, the frequency of supemovae, and the masses of 
compact remnants. The existence of a cutoff mass also has interesting implications for AY/AZ 
and the synthesis of ^^Al and ^°Fe which we also briefly discuss. 

We choose to analyze the statistics of the isotopes lighter and heavier than the iron group 
nuclei separately as we expect the elements heavier than atomic mass 60 to be secondary nucle- 
osynthesis due primarily to the 5-process or the reprocessing of 5-process isotopes in the original 
star. The mean and standard deviations of the production factors of the lighter and heavier iron 
group nuclei are tabulated separately in Table [T] 



3.1. The Production of Carbon, Oxygen, and Intermediate Mass Elements 



When studying the nucleosynthesis of massive stars, theorists often normalize to the produc- 
tion of ^^O, as it is the third most abundant isotope in the universe and comes almost entirely from 
massive stars (|Langerlll996|) . For reference, we've included the production factor of ^^O in Table [U 
as well. 

In Fig.[3]we compare the production factor of various elements to that of '^O for variable Mbh- 
For present purposes, success is defined as being within a factor of two of the solar abundances. 
We use a number of alpha elements to probe the strength of the oxygen and silicon burning, and 
we use ^"^N to measure the strength of the CNO cycle. From this, it's clear that little is lost in terms 
of the production of these common isotopes if the upper niass li mit is reduced from 120 to 40 Mq, 
which supports the preliminary limits set by iHeger et al.l (|2003|) . 

However, for low upper mass limits, '^C is overproduced in the winds compared to ^^O, 
suggesting either that this bound cannot be much lower than 25 M0 (a production factor of two) 
or that there is significantly less mass loss than assumed. Carbon is mostly made in the winds of 
very massive stars, especially during their Wolf-Rayet stage. Oxygen, which is the major part of 
the metallicity is also made in winds but more in the explosions. Their ratio then is quite sensitive 
to the mass loss prescription used. To illustrate this, in Fig. [3j we've also included the ratios for 
a population for which we've halved the mass loss yields. In the figure, it's apparent that the 
C/0 ratio remains consistent with the solar abundances down to an upper mass limit of 20 M0. 
To illustrate the change in the nucleosynthesis between the upper mass limit extrema, we present 
Fig- Has a comparison to Fig. |2] The apparent large overproduction of "^ "K in both figures is not 
a problem because "^"K is radioactive with a half life of 1.26 x 10^ years (|Lidell2009|) . Much of it 
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will decay before being incorporated into the Sun. 

We have checked other important intermediate mass isotope ratios in order to probe how 
well these stellar populations represent the solar abundances. These include ^*^Si/'*°Ca, '^°Ca/^'*Mg, 
28§j/56pg^ and ^°Ca/^^Fe. As can be seen in Fig. [51 these all remain acceptably close to solar ratios 
for the entire range of Mbh- 



3.2. The Supernova Rate 

Using the production factor of ^^O, we can estimate how the supernova rate depends upon 
Mbh- An accurate calculation of the rate itself would require a galactic chemical ev olution model 



that in cludes gas accretion onto and outflows from the galactic disk — as described in lTimmes et al 



(119951) — and is beyond the scope of this paper. However, we can estimate the factor by which the 
supernova rate would have to increase by examining Nsn, the total number of massive stars required 
to die to produce the correct amount of ^''O, normalized to the total number required to produce 
the yields of our control population (Mbh =120 M©). The number of stellar deaths (and therefore, 
supemovae) increases slowly with the lowering of the upper mass limit until approximately 40 
M0, below which it increases rapidly up to twice the original number at 28 M0 and three times the 
original number at 19 Mq. 



3.3. The Light 5^-Process 

The heavier component of the ^'-p rocess, those nuclei above A 90 is thought to be produced 
in low mass stars ( Pignatari et al.ll20rQ) . The light component of the i'-process, on the other hand, is 



generally attributed to massive stars and occurs late during helium burning when the temperature 
rises sufficiently for the ^^Ne(a,n)^^Mg reaction to occur. Owing to the temperature sensitivity 
of this rate, production of the i'-process isotopes in the mass range A = 60 - 90 is a potentially 
powerful constraint on Mbh- It is important that we oveq)roduce the 5-process elements somewhat 
in this stellar mass range as many of the primary isotopes are also made in lower metallicity stars 
that don't produce many 5-process elements. 

Use of this diagnostic is complicated by the fact that many isotopes in the mass range A = 60 
- 90 can be produced by both the s and r-processes, and so an alternative approach is to focus on a 
few "^•-only" is otopes (Fig. |9l)— thos e isotopes produced exclusively by the s-process — ™Ge, ^''Se, 



*^^Sr, and *^'Sr ([Kappeler et al.lll989|) . We exclude *^"Kr and ^^Kr from our analysis — despite both 
being 5-only isotopes — as these isotopes are also significantly produced in low-mass asymptotic 
giant branch stars. We notice that the 5-process is most predominant in core-collapse supemovae 
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from 30 to 50 M0. With more complex models of the stellar population, a comparison of the 
total 5-process yields to the solar abundance could prove to provide greater insight as to the upper 
mass bound. Here, we see that ^^Sr and ^^Sr are underproduced regardless of the upper bound, 
which could indicate an issue with nuclear cross sections or our stellar models. The other isotopes, 
^"Ge and ^^Se appear to be produced in sufficient quantities for us to use them for this paper. The 
production of ^°Ge falls below 1/2 the solar abundance at 21 M0, whereas the production of *^^Sr 
falls below 1/2 solar production at 23 Mfr^ . However, if we in crease the rate of ^^Ne(Q;,n) to the 
maximum experimental value according to IJaeger et al.l (|200l|) . we find that a limit of 18 Mq is a 
reasonable value for the mass of the heaviest supemovae. 



3.4. ^"Fe and ^^Al 



The nuclei ^"Fe (ri/2 = 2.62 x 10^ y) and ^^Al (ri/2 = 7. 17 x 10^ y) are interesting because they 
accumulate in the interstellar medium where t he emis s ion ge nerated by their decays can be studied 
using gamma-ray telescopes. Observations by ISmithI (|2004|) give a ratio for the decay rate of ''^Fe 
to that of Al of about 1 6 implying a ratio of ^"Feto^^Al m ass fractions of (6 0/26)(0.16) = 0.37 
( Wooslev & Hegeill2007h . As pointed out by IPrantzosI (|2004() and discussed in IWoosley & Heger 
(120071) . current calculations give a larger value. Our c urren t study uses yields that already address 
some of the concems discussed in lWoosley & Hegeii(|2007|) . including inappropriate rates for ^^Al 
destruction and use of OPAL opacities for electron scattering at high temperature. As Fig.|6]shows, 
our mass-averaged estimate of the production is 1.0 for Mbh = 120. This compares favorably with 
the value 0.95 given for these models by the authors. 

The remaining difference between 1.0 and 0.37 probably chiefly reflects remaining uncer- 
tainty in the neutron capture cross sections for ^^Fe and ^^Al and especially the choice of ex- 
plosion energies for low mass supemovae. Studies of the light curves of Type Hp supemovae 
(|Kasen & Woosleyll2009l) . suggest that many supemovae, perhaps most, have an explosion energy 
smaller than the fiducial 1.2 x 10^' erg assumed in many st u dies. Reducing the explosion energy 
reduces the yield of ^°Fe substantially. iLimongi & Chieffil (|2006|) have also emphasized the de- 
pendence of this production ratio on mass loss, convection theor y, and the initial ma ss function. 
Including the effects of rotation may also increase the ^^Al yield (iPalacios et al.ll2005|) . especially 
in very massive stars. 

Our purpose here is not to completely solve the debate surrounding ^°Fe/^^Al but to point 
out that the answer depends upon Mbh- Measurements of gamma-ray line flux ratios may thus 
ultimately help constrain the masses of stars that explode. Fig. [6] shows that the ratio of ^"Fe/^^Al 
produced by supemovae declines as Mbh becomes smaller. Given that the production of ^"Fe will 
be even smaller in those stars we are assuming still explode if their kinetic energy is reduced, this 
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effect could reduce the average production of Fe appreciably. 



3.5. Helium Production and AY/AZ 



As discussed extensively in lMaeden (|1992|) . the measured derivative of the helium mass frac- 
tion with respect to metallicity can, in principle, be used to constrain Mbh- This is because the 
winds of the most massive stars are rich in helium while the heavy elements are largely confined 
to the cores. If the cores collapse to black holes, trapping the heavy elements, the synthesis of the 
winds remains and increases the overall average of AY/AZ. Observations suggest a value AY/AZ 
~ 4. 

In practice, the application of this metric is fraught with uncertainty. The yields of the massive 
stars are dependent upon the initial mass function and especially the very uncertain mass loss rates 
employed. After the helium core is uncovered by mass loss, Wolf-Rayet mass loss contributes not 
only helium, but an increasing amount of carbon and oxygen. The remnant masses depend upon 
an uncertain explosion mechanism. Does all the presupemova star fall into the hole or only part? 
For successful explosions, where is the mass cut? Lower mass stars also contribute appreciably 
to both helium and metallicity and the yields of stars of all masses are sensitive to metallicity. 
Even an approximate meaningful result requires integration over some uncertain model for galactic 
chemical evolution. 

However, because this metric has been applied extensively in the literature, we give in Fig. |7] 
our results for the solar metallicity stars considered in this survey. Until Mbh is reduced below 
~40 M0 AY/AZ remains slightly less than unity. Even for Mbh = 18 Mq, AY/AZ is still only 
1.929, well below the observed value. 



These results are, at first glance, seemingly inconsistent with those of Timmes et al.l (|1995|) 



who found AY/AZ = 4 for Mbh = 17 Mq for sola r metallicity stars and 30 M f;^ for low metallicity 



stars. However, Timmes et al. used the survey of IWoosley & Weaver! (| 1995b which included only 



stars below 40 M©, while our grid extends to 120 Mq. Furthermore, Timmes et al. assumed that 
any winds would be metal free and not extend into the helium core. We include initial metals 
in the envelope and mass loss from Wolf-Rayet stars once the core is uncovered. If we redo our 
calculation using the assumptions of Timmes et al., the dashed curve in Fig. |7]results. This curve 
crosses AY/AZ at Mbh = 17 Mq, in excellent agreement with Timmes et al. However, we believe 
that our results for solar metallicity are more realistic and that a lower value of AY/AZ from 
massive stars is appropriate. 
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3.6. Average Mass of Compact Remnants 



We use Table 4 from IZhang et al.l (|2008l) to get the baryonic masses of the remnants for our 
fiducial set of supemovae from stars of initial masses of 10 to 100 Mq. We then calculate the 
gravitational mass of the stellar remnants with their Equation 2. We assume that the maximum 
gravitational mass of a neutron star is 2.0 Mq so all remnants with a higher gravitational mass 
than this form black holes with a gravitational mass equal to its baryonic mass. We then construct 
a population of stars according to a Salpeter initial mass function and vary the heaviest supernova 
mass. We let any star with M > Mbh become a black hole with mass equal to its presupemova 
mass. In Fig. [8l we plot various quantities (the average remnant mass, the average black hole 
mass, the average neutron star mass, the average iron core mass, and the maximum remnant mass) 
as a function of Mbh- For those stars more massive than the heaviest supernova, we have assumed 
that the presu pemova star compl etelycoUapses to a black hole, so our black hole masses are larger 
than those of Izhang et al.l hoO^ . 



We can compare this to the observed average neutron star mass from ISchwab et al.l (|2010|) . 
From 14 neutron stars with well-measured masses, they find an average neutron star mass of 
1.325 ±0.056 Mq. We find an average neutron mass closest to this at 14 Mq and are within 
their uncertainty up to 16 Mq. This suggests that a low value for Mbh is more consistent with the 
observed neutron star masses. 



4. Conclusions 

By examining a variety of nucleosynthetic diagnostics, one can constrain the mass of the max- 
imum mass supernova that must explode and not swallow up its heavy el ements in a black hole. A t 



the outset, a number of caveats are worth stating. First, our yield set from lWoosley & Hegeii (|2007|) 
reflects a particular choice of many uncertain parameters — the treatment of convection, mass loss, 
key reaction rates (like ^^C(q!,7)'^0), and explosion physics. We lump together both means of 
forming a black hole, direct and by fall back and we ignore mixing. It is quite possible, probable 
even, that some supemovae make black holes and yet still eject some fraction of their presuper- 
nova nucleosynthesis. Mixing during the fallback epoch can complicate the estimate of yields. 
This probably affects the intermediate mass elements and iron more than helium, carbon, oxygen 
and the i'-process. Our models sample a single metallicity — solar. While other models of different 
metallicity are available, using them would require a more careful treatment of galactic chemical 
evolution than is presently justified, given all the other uncertainties. Metallicity will not greatly 
affect the supemova yield of primary elements like oxygen and the intermediate mass elements, 
but it does affect the mass loss. We attempted to partially compensate for that by multiplying the 
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integrated nucleosynthesis of the wind by a constant factor of 0.5, as shown in Fig. [31 The use 
of solar metallicity also overestimates somewhat the pro duction of secondar y species such as the 



elements with odd nuclear charge and the i'-process (e.g. jTimmes et al.lll995h . so somewhat larger 



yields than solar might be required than suggested by Fig. |4] Finally, all our models ignore rota- 
tion. Not only can rotation affect the mechanism and symmetry of the explosion, but it also makes 
the mass of the helium and carbon-oxygen core larger for a given main sequence. So our derived 
mass limits may actually need to be smaller. 

Given these limitations, the best we can say at the present time is what supernova mass limits 
might be consistent with observations. The idea of a limiting mass is itself an appr oximation, since 



the co mpactness of the core is not a monotonic function of main sequence mass ( O'Connor & Ott 



2011|) . especially in the interesting range 20 - 35 Mq. For still heavier stars, mass loss may shrink 
the helium core so much that the presupemova helium core mass of say a 100 Mq star differs little 
from that of a 20 M© star. Such massive stars are rare however, and their nucleosynthesis is mostly 
due to their presupemova winds. 

We have looked at several processes that limit Mbh- As Mbh is decreased, the necessary 
rate of "successful" supemovae rises. For Mbh = 20 the rate is 2.88 times greater than for Mbh 
= 120. Surprisingly, in reducing Mbh to 20 Mq, the overall nucleosynthesis of most isotopes 
from Ne to Ca with respect to '^O is altered little and the production of some isotopes in the iron 
group is actually improved. The greatest apparent problem is ^^Ne; however, this can be mitigated 
somewhat by reducing the amount of mass loss or by increasing the ^^Ne(a,n) rate. Of the s- 
process only isotopes produced mainly in massive stars, ™Ge becomes the limiting isotope for the 
mass of the heaviest supernova, reaching half of solar value at 23 Mq. This can be extended to 18 
Mq by increasing the ^^Ne(a,n) rate. 

In total, we find that we can easily reduce the mass of the heaviest supemovae to 40 Mq 
without any significant changes. For KEPLER's standard values of nuclear rates and mass loss, 
there are only moderate changes in these processes down to 25 Mq. The limits become increasingly 
severe for smaller masses. For Mbh = 25 Mq, the stellar winds overproduce ^^C with respect to 
^^O by a factor of two, unless we reduce the mass loss in these stars by two. At Mbh = 21 Mq, 
the lighter elements are overproduced in these more massive stars, and the 5-process produces only 
half the needed ^"Ge unless the ^^Ne(Q;,n)^^Mg rate is increased. At 20 Mq, even with halved mass 
loss, the winds overproduce ^^C by a factor of two. 
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Fig. 1. — The Salpeter Initial Mass Function is a model that represents the distribution of stellar 
mass in a population. Mass here is plotted logarithmically. The binning indicates the binning used 
in our grid of stellar models. 
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Fig. 2. — The production factor (plotted logarithmically) as a function of atomic mass. The dashed 
line represents the production factor of '^O. The dashed lines represent a factor of two deviation 
from ^^O. 
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Fig. 3. — The relative abundances of the common isotopes with respect to ^''O as a function of Mbh- 
The plot has been truncated at Mbh = 6OM0, beyond which the abundances do not illustrate any 
significant variations. A relative abundance of 1 indicates solar ratio. The dotted lines represent 
the ratios produced if the wind contribution is halved to illustrate the dependence on wind. 
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Isotopic Abundances (Upper Mass Limit 20) 
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Fig. 4. — The production factor (plotted logarithmically) as a function of atomic mass. The dashed 
line represents the production factor of ^^O. The dotted lines represent a factor of two deviation 
from ^^O. Here it is assumed that all massive star over 20 M0 end up collapsing as black holes 
and their only nucleo synthetic contribution is their winds. The large abundance of "^^K is not a 
problem as this species will decay prior to solar system formation bu the large abundances of ^^Ne 
and ^^C could indicate unresolved issues with mass loss, nuclear physics, or the assumption of a 
low threshold for black hole formation. The ^-process is also weak here (compare with Fig. |2]). 
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Fig. 5. — The relative abundances of the intermediate mass isotopes as a function of Mbh- The plot 
has been truncated at Mbh = 60Mq, beyond which the abundances do not illustrate any significant 
variations. 
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Fig. 6. — The ratio of the total masses of ^'^Fe to ^^Al as a function of upper mass Umit. 
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Chemical Enrichment 
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Fig. 7. — The ratio of the total deviations from the initial helium and metal mass fractions of the 
population as a function of the upp er mass limit. The curve labeled "Truncated" uses the same 
assumptions as lTimmes et al.l (| 19951) in order to recreate the solar metallicity curve of Figure 37 of 
that same paper. 
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Average Remnant Mass 
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Fig. 8. — The average neutron star remnant mass, iron core mass, average black hole mass, max- 
imum remnant mass, and average remnant mass are plotted as a function of Mbh- The average 
neutron star mass is within the observational uncertainty everywhere below Mbh = 25 M0. 
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Relative Abundances of S-Only Isotopes 
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Fig. 9. — The relative abundances of the s-only isotopes with respect to °0 and "^Mg as a function 
of Mupper- The plot has been truncated at Mupper = 6OM0, beyond which the abundances do not 
illustrate any significant variations. A relative abundance of 1 indicates solar ratio. The dashed 
lines represent the ratios produced if the ^^Ne(Q!,n) is increased to the maximum experimental 
value. 
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Table 1. Production Factor Statistics 



Mupper (M©) 


l6r» 
(J 


Mean (io-jy) 


a (io-jy) 


Mean (oU-o4) 


a (,oU-o4) 




120 


14.74 


6.21650 


5.48443 


14.59968 


6.72117 


1.0000000 


100 


14.63 


6.18653 


5.45782 


14.54773 


6.69476 


1.0149254 


80 


14.48 


6.14438 


5.42532 


14.48140 


6.66672 


1.0298507 


70 


14.24 


6.04946 


5.37006 


14.34487 


6.61504 


1.0447761 


60 


14.01 


5.96336 


5.30620 


14.16470 


6.54938 


1.0597014 


50 


13.14 


5.60242 


4.78830 


13.20545 


6.08493 


1.1343282 


40 


11.46 


5.19207 


4.36852 


10.63252 


4.68389 


1.2985075 


30 


7.93 


4.23886 


3.51987 


5.86116 


2.37390 


1.8805970 


28 


7.33 


3.97234 


3.26334 


5.26788 


2.19142 


2.0298507 


25 


6.52 


3.60184 


2.96335 


4.35242 


1.87591 


2.2835820 


24 


6.27 


3.49595 


2.91884 


3.95001 


1.78152 


2.3731341 


23 


5.99 


3.37648 


2.87740 


3.59129 


1.71603 


2.4925373 


22 


5.69 


3.25718 


2.84823 


3.29207 


1.60152 


2.6119401 


21 


5.44 


3.14271 


2.82677 


3.01019 


1.51820 


2.7462685 


20 


5.17 


3.00155 


2.81088 


2.69667 


1.46384 


2.8805971 


19 


4.93 


2.82092 


2.76516 


2.51345 


1.39076 


3.0149255 


18 


4.70 


2.66811 


2.78619 


2.30015 


1.31519 


3.1791043 


17 


4.52 


2.55969 


2.81806 


2.11734 


1.16281 


3.2985075 


16 


4.36 


2.43725 


2.87535 


1.94110 


1.00387 


3.4179103 


15 


4.22 


2.31317 


2.96168 


1.78076 


0.83437 


3.5223880 


14 


4.12 


2.20852 


3.07487 


1.60545 


0.62216 


3.6119401 



